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ABSTRACT 

It is widely accepted that feedback from active galactic nuclei (AGN) plays a key role in the evolution 
of gas in groups and clusters of galaxies. Unequivocal evidence comes from quasi-spherical X-ray 
cavities observed near cluster centers having sizes ranging from a few to tens of kpc, some containing 
radio emission. Cavities apparently evolve from the interaction of AGN jets with the intracluster 
medium (ICM). However, in numerical simulations it has been difficult to create such fat cavities from 
narrow jets. Ultra-hot thermal jets dominated by kinetic energy typically penetrate deep into the 
ICM, forming radially elongated cavities at large radii unlike those observed. Here, we study very 
light jets dominated energetically by relativistic cosmic rays (CRs) with axisymmetric hydrodynamic 
simulations, investigating the jet evolution both when they are active and when they are later turned 
off. We find that, when the thermal gas density in a CR-dominated jet is sufficiently low, the jet has a 
correspondingly low inertia, and thus decelerates quickly in the ICM. Furthermore, CR pressure causes 
the jet to expand laterally, encounter and displace more decelerating ICM gas, naturally producing 
fat cavities near cluster centers similar to those observed. Our calculations of cavity formation imply 
that AGN jets responsible for creating fat X-ray cavities (radio bubbles) are very light, and dominated 
by CRs. This scenario is consistent with radio observations of Fanar off- Riley I jets that appear to 
decelerate rapidly, produce strong synchrotron emission and expand typically at distances of a few 
kpc from the central AGN. 

Subject headings: cosmic rays - galaxies: active - galaxies: clusters: intracluster medium - galax- 
ies: jets - X-rays: galaxies: clusters 



1. INTRODUCTION 

Feedback from active galactic nuclei (AGNs) is increas- 
ingly recognized to play an important role in the evolu- 
tion of elliptical galaxies, groups and clusters of galax- 
ies. It is considered to be a prime candidate for solv- 
ing the "overcooling" problem associated with galax- 
ies at the brig ht end of the galaxy luminosity function 
(jCroton et all 120061 : iTeyssier et all I2Q1QD . i.e., cosmo- 
logical simulations without AGN feedback predict mas- 
sive galaxies that are too bright and too blue when 
compared to massive galaxies in the nearby u niverse 
(|Kravtsov et all 120051 : feorgani fc Kravtsovll2009f ). Mas- 
sive galaxies sit at the centers of large dark matter haloes, 
and the "overcooling" problem is thus related to the so- 
called "cooling flow" problem in galaxy clusters, groups, 
and giant elliptical galaxies. High-resolution X-ray ob- 
servations from the Chandra and XMM-Newton Tele- 
scopes have shown that, though the gas in these systems 
emits copiously in X-rays, the amount of gas cooling to 
very low temperatures is much less than the amount of 
hot g as with a short cooling time (see lPeterson fc Fabianl 
2006 for a recent review). This suggests that some heat- 
ing mechanism is at work to counterbalance radiative 
losses. It is now widely agreed that AGN feedback plays 
a key role in heating the h ot gas in these systems (see 
McNam ara fc Nu lsen 200.3 for a recent review). 

The strongest observational evidence for AGN feed- 
back is in galaxy clusters where high-resolution X- 
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ray observations have detected numerous X-ray-deficient 
cavities, some of which are associated with radio 
j ets and spatially coinciden t with radio bubbles/lobes 
([McNamara fc Nu lsen 2007). Most cavities have spatial 
sizes of a few to a few tens kpc and are located near clus- 
ter center s with cluster-centric distan c es of a few to a few 
tens kpc (jMcNamara fc Nulsenll2007l : iDiehl et al.ll2008D . 
These cavities are roughly spherical or elongated in the 
tangential direction (i.e., perpendicular to the radial jet 
axis), and are likely rising buoyantly in the intraclus- 
ter medium (ICM). Observations indicate that roughly 
70 - 75% of X-r ay bright cool-core clusters harbor de- 
tectable cavities (jDunn et al.ll2005l ), and the associated 
AGN energies are sufficient or nearly sufficient to stop the 
gas from cooling in most clu sters containing detectable 
cavities ([Raffertv et all 12006( 1. This suggests that AGN 
feedback indeed plays a significant ro le in the thermo - 
dynamical evolution of galaxy clusters (Guo et al. 2008). 
X-r ay cavities have al so been detected in elli ptical galax- 
ies (| Allen et al.ll2006[ ) and in galaxy groups ([Dong et al.l 
1201 Of ), but it is usually more difficult to detect cavities 
in smaller halos with fainter X-ray emission. 

Combined radio and X-ra y observations of X-ray cavi- 
ties (e.g., IDunn et al.l 120051 : [McNamara et all 12009) sug- 
gest that these cavities are inflated by bipolar jets ema- 
nating from AGNs in cluster-central galaxies. It is natu- 
ral to study AGN feedback by investigating the evolution 
of AGN jets in these systems. In recent years, numerical 
simulations of jets in gal axy clusters have be e n conducted 
by many authors (e.g., [O mma et al. 2004; Zanni et al.l 
120051: IVer naleo & Revnolds 2006; BruggeneFaD 120071 : 
O'Neill & Jones 2010; Morsony et al. 2010), who mainly 
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focus on the impact of AGN outbursts on the thermal 
and dynamical evolution of the ICM. In contrast, there 
are relatively few gas-dynamical studies of the physics 
of AGN outbursts itself or the content (composition) of 
AGN jets that directly affect the impact of AGN feed- 
back on the ICM. In most of these calculations, it is 
assumed that AGN jets contain ultra-hot thermal gas 
with jet-to-ICM density contrasts r] = Pj/p ~ 0.01. 
The energy fluxes of these thermal jets are commonly 
dominated by the kinetic energy. Simulations indicate 
that these jets usually penetrate deep into the ICM, 
forming radiall y-elongated cavities at large radii (e.g., 
Reynolds et~aTl I2002L lOmma et all 12004 IZanni et all 
2005L iVernaleo fc Reynold s 2006. O'Neill fc Jonesl l2010l) . 
but never forming the "fatter" quasi-spherical X-ray cav- 
ities observed near cluster centers. This 'cavity forma- 
tion' problem suggests that the adopted model for either 
AGN jets or the ICM has at least one inaccurate key 
ingredient, which may lead to quantitively or even qual- 
itatively incorrect results. 

The c avity formation problem has been previously no- 
ticed bv lSternberg et aTl (|2007[ ). who propose wide, con- 
ical jets with typical half-opening angles a > 50° to 
form fat cavities in galaxy clusters. Rapidly precess- 
ing jets with large precessin g angles can also produce 
quasi-spherical fat cavities ([Sternberg fc Sokerf [2008a; 
but slowly precessing jets may produce multiple pairs 
of cavities, Falceta-Goncalves et al. 2010). Alternatively, 
jet simulations by Briigge n et al.l ((2007) suggest that nar- 
row jets could result in fat cavities if significant random 
motions are pres ent in the ICM . Fat c avities are seen 
in simulations of iMorsonv et al.l (j2Q10h , which include 
both jet precession and ICM random motions. How- 
ever, it is unclear if any of the above three conditions is 
common in extragalactic AGN jets responsible for form- 
ing fat cavities. Consider for example the Virgo cluster, 
the closest galaxy cluster which exhibits de tailed interac- 
tions of AGN feedback with the ICM (e.g., Mil lion et al. 
2010). The 90 cm radio image of M87 bv 10 wen et al. 
(2000]) shows two regular, nearly spherical radio lobes, 
unlikely to have been strongly shaped by random gas 
motions which tend to form ir regular bubbles (see Fig- 
ure 1 in Briigge n et al.l 120071 ). Multi- wavelength ob- 
servations of the famous M87 jet clearly s how that it 
is highly collimated even on parsec scales ([Junor et al.l 
1999; iKovalev et al.l 12007). Narrow radio jets have also 
been observed in numerous other systems, indicating 
that the jet opening angle is generally very small. If the 
eastern ear-shaped radio lobe is formed by t he jet traced 
by th e eastern radio arm (see Figure 2 in 10 wen et al.l 
2000), the well-collimated eastern arm may suggest that 
th e jet does not have a large p recession angle as required 
bv lSternberg fc Sokerl (|2008al ) to create fat bubbles. Ob- 
servations of M87 thus motivate us to explore other miss- 
ing ingredients of AGN feedback, which may be respon- 
sible for forming fat cavities. 

There is indeed one important ingredient of AGN feed- 
back missing in previous simulations - the relativistic 
cosmic ray (CR) component of AGN jets. CRs provide 
pressure with negligible inertia, which can puff up jets 
to form distended cavities and radio lobes near cluster 
centers. The same effect can also be achieved by thermal 
jets with very low density contrasts and high internal 
pressures, but thermal particles in these jets are actually 



relativistically hot (see discussions in Section lcT3|) . Fur- 
thermore, numerous detections of radio jets and lobes as- 
sociated with X-ray cavities point toward the existence of 
a significant CR electron population. Deep radio obser- 
vations of Fanar off- Riley type I (FR I) radio jets, which 
are likely to be responsible for creating X-ray cavities, 
indicate that these jets decelerate rapidly and produce 
strong synchrotron emission (flaring) a t typical distances 
of a few kpc from central nuclei (e.g., lLaing et al.l 2006; 
lLaing et al.l 12008). For the first time, in this paper we 
perform numerical simulations of CR-dominated jets and 
follow the dynamical evolution of these jets in a typical 
galaxy cluster. These jets contain a dominant CR com- 
ponent and a very small thermal component. They are 
very light and yet have a large energy flux and high in- 
ternal pressure. Due to the reduced jet inertia and mo- 
mentum with enhanced lateral expansion driven by CR 
pressure, these jets naturally produce fat cavities near 
cluster centers, resembling observed cavities very well. 

Radio jets have been previously proposed to explain 
extended (part i cularl y FR II) e xtragala ctic radio sources 
(lLongair et al.l 119731 : IScheuerl 119741 : iBlandford fc Reesl 
119741 ). The propagation of purely relativistic radio 
j ets has been studied with conserva tion-law calculations 
(jBicknelllll99l lLaing fc Bridldl2002D . but numerical sim- 
ulations usually focus on jets composed of ideal thermal 
gas (some incl uding magnetic fields o r relativistic bulk 
motions, e.g., Perucho & Marti 2007, O'Neill & Jones 
I2010D . Unlike IScheck et all (|2002f ) and iPerucho fc Martil 
(2007), who adopt an equation of state for a thermal- 
ized relativistic ideal gas, our hydrodynamical calcula- 
tions follow an additional relativistic CR component that 
can have a power-law energy distribution as expected 
from shock acceleration. We find that fat X-ray cavities 
observed in galaxy groups and clusters are most natu- 
rally formed from jets in which the CR energy density 
dominates and which are very light, i.e. containing non- 
relativistic thermal gas with densities that are ~ 10 -4 
times the central cluster gas density. Such light jets have 
also been con s idered for example by iKrausd ([2003) and 
iGaibler et al.l (|2009[ ) in studies of purely thermal FR II 
jets and cocoons. For the first time, we numerically in- 
vestigate the evolution of two-fluid jets composed of both 
thermal gas and CRs. 

The rest of the paper is organized as follows. In Sec- 
tion [2j we describe the physics of galaxy clusters with 
CRs, and numerical setup. We present our results, in 
particular, how CR-dominated jets form fat cavities, in 
Section [3j and summarize our main conclusions with im- 
plications in Section HI 



2. EQUATIONS AND NUMERICAL METHODS 

2.1. Cluster Physics with Cosmic Rays 

X-ray cavities are usually thought to be inflated by 
bipolar jets emanating from AGNs located at cluster cen- 
ters. Radio synchrotron emission has been detected from 
many jets and X-ray cavities, confirming the existence of 
a significant CR component, which may be transported 
along the jets and/or created in strong shocks as the jets 
encounter the ICM. The combined hydrodynamic evo- 
lution of thermal gas and CRs can be described by the 
following four equations: 
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Fig. 1. — Central slices (100 x 20 kpc) of log (n e /cm 3 ) in runs JO (top panels), Jl (middle panels) and Jl-A (bottom panels) at time 
t = 10 Myr (left panels) and t = 70 Myr (right panels). The horizontal and vertical axes represent z and r respectively. 
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where d/d£ = d/dt + v • V is the Lagrangian time deriva- 
tive, P c is the CR pressure, e c is the CR energy density, k 
is the CR diffusion coefficient, S c is the CR source term, 
and all other variables have their usual meanings. Pres- 
sures and energy densities are related via P = (7 — l)e 
and P c = (7 C — l)e c , where we assume 7 = 5/3 and 
7 C =4/3. We ignore radiative cooling, which is unimpor- 
tant during our short-duration (< 100 Myr) simulations 
of cavity formation. The gas temperature is related to 
the gas pressure and density via the ideal gas law: 



k B p 



(5) 



where &b is Boltzmann's constant, m M is the atomic mass 
unit, and \i = 0.61 is the molecular weight. To avoid 
confusion we denote the gas pressure P as P g in the rest 
of the paper. 

Our model and methods are generally applicable to 
all reasonably relaxed clusters, but for concreteness, we 
adopt simulation parameters appropriate for the well- 
observed nearby Virgo cluster. The initial radial profiles 
of gas density and temperature are take n from the ana- 
lytic fits to the observations provided by Ghizzardi et al. 
(|2004f ). We choose a gravitational field that establishes 
exact hydrostatic equilibrium for the initial gas pressure 
profile. At time t = we assume that the CR energy den- 
sity is zero in the cluster gas, e cr = 0, but at later times 
CRs enter the cluster in jets. Equation 2] describes the 
CR advection and diffusion in the ICM. The CR diffusion 
coefficient n is poorly known but may vary inversely with 



the gas densit y since the magnet ic field is probably larger 
in denser gas (Dol ag et al.ll200lf ). We adopt the following 
functional dependence of the diffusion coefficient on the 
electron number density n e : 



10 29 (n e o/^ e ) cm 2 
10 29 cm 2 s- 1 



when n e 
when n e 



> ™e0, 
< ™e0, 



(6) 



where n e o = 10 -5 cm -3 . This level of diffusion does 
not strongly affect the early evolution of cavities, and 
our results are fairly insensitive to it. During their dif- 
fusion, CRs interact with magnetic irregularities and 
Alfven waves, effectively exerting CR pressure gradi- 
ents on the thermal gas (equation [2j). Pressure gradi- 
ents in the CR component act directly on the gas by 
means of magnetic fields frozen into the gas but the 
magnetic energy densities are too small to contribute to 
the gas dynamics and do not appear in the equations 
above. Recent observations also suggest that the mag- 
netic field is likely to be dynamically subdominant in 
X-ray cavities, where the energy density is mainly con- 
tributed by a particle population that is not responsi- 
ble for the observed synchrotron emission ([Croston et al.l 
l2CMlHaxdcastle fc CrostoiJ2QlrtlDunn et al.ll2010l ). We 
neglect other more complicated (probably secondary) in- 
teractions of CRs with thermal gas, e.g., Coulomb in- 
teractions, hadronic collisions, and hydromagnetic- wave- 
mediated CR heating, that depend on the CR energy 
spectrum and may provide additional heating effects for 
the ICM (e.g., IGuo fc Ohll2008h . We also neglect CR 
energy losses from synchrotron emission; although syn- 
chrotron losses may be important for high-energy CR 
electrons, most of the contribution to the integrated CR 
energy density e c comes from low-energy CR electrons 
and possibly CR protons. 

Our method of CR modeling has been successfully 
used in a series of papers, e.g., [ Mathews & Brig hentH 
(2008a), Mathews & Brighenti (2008b), Mathews (20C_ 
IGuo fc M athews (2010b), and IGuo fc Mat hews (2010a 
to which the reader is referred to for further details. 
When X-ray cavities were formed in these previous stud- 
ies, we ignored the complex jet physics and assumed that 
CRs were injected in situ into a small region slightly off- 
set from the cluster center. Once formed, X-ray cavities 
rise buoyantly in the ICM. In this paper, we take a step 
further and study in detail how CR-dominated jets form 
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cavities in the ICM. To this end, we assume that the CR 
source term in Equation 3] is S c = 0, and all the CRs are 
injected into the ICM by the jet at the inner boundary of 
our simulations, as explained in the following subsection. 

2.2. Numerical Setup and Jet Injection 

Equations (pp) — (|4|) are solved in (r, z) cylindrical co- 
ordinates usin g a two-dimensional Eul erian code similar 
to ZEUS 2D ([Stone fc Normanl IT992). Our code con- 
tains a background gravitational potential, CR diffusion, 
and an energy equation for CRs. The computational 
grid consists of 400 equally spaced zones in both coordi- 
nates out to 100 kpc plus additional 100 logarithmically- 
spaced zones out to 1 Mpc. For both thermal and CR 
fluids, we adopt outflow boundary conditions at the outer 
boundary and reflective boundary conditions at the inner 
boundary except for those inner boundaries representing 
the jet base when the jet is active. The jet inflow is in- 
troduced along the z-axis through a nozzle placed at the 
cluster center. At its base, the jet has a radius of rj e t = 1 
kpc (4 zones along the r axis) and is initi alized by inflow 
bound ary conditions similar to those in Briig gen et al.l 
(|2QQ7h . By preparing the jet properties before it enters 
the computational grid, we avoid an initial entrainment 
of cluster gas into the jet as usually occurs if the jet is 
initialized within active zones. 

The jet contains both thermal gas and CRs and is in- 
jected with cylindrical geometry having an initial open- 
ing angle of degrees. At their base in central ghost 
zones jet inflows are described by six parameters: the 
density ratio n = p^/po of thermal gas within the jet to 
the ambient gas at the cluster center, the energy density 
of thermal gas within the jet e^, the CR energy density 
within the jet ej cr , the jet velocity Uj et , the jet radius 
rjet, and the jet duration tj e t- The jet power can then be 
written as 

Pjet=Pke+Pcr + Pth, (7) 

where Pk e = pjV? et irr? et /2 is the jet kinetic power, P cr = 
e jcr^jet7rr?et is the jet CR power, and P th = e^v^ et 7rrf et 
is the jet thermal power. The total energy injected by 
the jet can be written as Ej e t — Pjet^jet- In the com- 
putations presented in this paper, we take Vj et = 0.05c, 
where c is the speed of light, rj e t = 1 kpc, and tj e t = 10 
Myr (except that in run J3, tj et = 2.72 Myr to keep 
Ejet the same as in our main run Jl). The jet has a 
Mach number of ~ 23 with respect to the initial cluster 
gas at the origin. AGN jets on parsec scales likely have 
relativistic velocities, but they may und ergo rapid decel- 
eration on kpc scales (jLaing et al.ll2006h . Our initial jet 
speed Vj e t = 0.05c represents jets that have already gone 
through the initial deceleration stage. We have also run 
simulations with Vj et = 0.1c or 0.2c, and found quali- 
tatively similar results. The jet internal thermal energy 
density ej is determined by n and the jet thermal temper- 
ature Tj, which is chosen to be 100 times the initial ambi- 
ent gas temperature at the origin. Since the jet power in 
our simulations is dominated by either the kinetic power 
or the CR power, our results are not sensitive to the value 
of Tj unless it is so high that ej > e- ]CT . The jet power in 
our main run Jl is Pj e t ~ 1.3 x 10 44 erg s _1 , which corre- 
sponds to a powerful FR I radio source and is chosen to 
roughly produce X-ray cavities with radii ~ 10 kpc near 



TABLE 1 
List of Simulations 



Run 


V 


cicr 


E C Y b 


E ke c 


E th d 


^jet 6 


*V 


JO 


0.01 








3.13 


1.07 


4.21 


1.3 


Jl 


0.0001 


2.96 


4.17 


0.03 


0.01 


4.21 


1.3 


Jl-A 


0.0001 








0.03 


0.01 


0.04 


0.01 


J2 


0.0001 


1.52 


2.14 


0.03 


0.01 


2.18 


0.69 


J39 


0.01 


8.0 


3.06 


0.85 


0.29 


4.21 


4.9 



a The initial CR energy density at the jet base (in units of 
10 -9 erg cm -3 ). 

b E C r, E ke , E t h, and Ej e t are, respectively, the injected CR, 
kinetic, thermal, and total energy (in units of 10 58 erg) by the 
jet during its active phase t < tj e t ( e -g- 5 E cr = P cr £j et ). 

•^The jet power in units of 10 44 erg s -1 . 

9 ln run J3, the jet is active for a duration of tj et = 2.72 Myr 
(to keep .Ejet the same as in run Jl), while in all other runs, 
tjet = 10 Myr. 

the center of our fiducial cluster. The jet energy and 
power are consistent with typical values es timated from 
cavit y energetics, which span a large range ([Birzan et al.l 
2004). In observation, the AGN mechanical luminosity is 
usually averaged within the cavity rising time (the buoy- 
ancy timescale), which is likely much longer than the jet 
duration. See Table Q] for the jet parameters and ener- 
getics in all runs presented in this paper. 

3. RESULTS 
3.1. Kinetic Energy- dominated Jets 

Numerical simulations suggest that jets dominated by 
kinetic energy usually penetrate easily into the ICM, 
forming low-density cavities at cluster centric distances 
much larger than those in real cluster s observed at 
similar ages (e.g., [Revnolds et al. 2002, lOmma et al.l 
2004. iZanni et all I2005L iVernaleo fc Reynolds! I2006L 
O'Neill fc Jones! 120101 ). Furthermore, these jet-inflated 
cavities are strongly elongated in the (radial) jet direc- 
tion, while most observed cavities are quite spherical, or 
elongated in the tange ntial direction (e.g ., the cavities 
in the Perseus cluster, iFabian et al.l 120061 ). The goal of 
this paper is to compute the dynamical evolution of CR- 
dominated jets and investigate if these physically moti- 
vated jets can form fat X-ray cavities near cluster centers. 

We first present our control run JO, which follows the 
evolution of a typical kinetic energy-dominated jet. As 
listed in Table [TJ the parameters of run JO are n = 0.01, 
and ej cr = 0, similar to those often adopted by previous 
authors (e.g., Revnolds et al. 2002, Vernale o fc Revnoldsl 
120061 . lO'Neill fc Jones! I201Q[ ). As usual, we assume that 
the jet thermal temperature at the jet base is 100 times 
the initial cluster gas temperature at the origin. Thus the 
jet is initially in pressure equilibrium with the ambient 
gas. This assumption is the most used in jet simulations, 
but may not hold for supersonic jets moving down the 
cluster pressure gradient. Furthermore, the initial jet 
pressure is determined by physical processes directly as- 
sociated with the central supermassive black hole, and is 
unlikely to be equal to the pressure of central hot ICM. 
The shock-heating of the jet material and the production 
of CRs within the jet may make the jet over-pressured. 
We will explicitly explore the effect of jet pressure on the 
formation of fat cavities in Subsection 13.31 

Figure [T] shows central 2D slices (100 x 20 kpc) of elec- 



COSMIC RAY-DOMINATED AGN JETS 



5 




Fig. 2. Same as FigureQ] but at t = 100 Myr. 

tron number density in logarithmic scale in runs JO (top 
panels), Jl (middle panels) and Jl-A (bottom panels) at 
time t = 10 Myr (left panels) and t = 70 Myr (right pan- 
els). The same images at t = 100 Myr are shown in Fig- 
ure O As clearly seen in these figures, the kinetic-energy- 
dominated jet in run JO penetrates easily into the ICM, 
and reaches ~ 100 kpc at t = 100 Myr. The jet creates a 
low-density cavity, which is strongly elongated in the jet 
direction. These features have been previously seen in 
many simulations of kinetic-energy-dominated jets, and 
are clearly inconsistent with observations of most X-ray 
cavities, which are often approximately spherical and ap- 
pear to rise buoyantly at cluster centric radii < 50 kpc 
in times t < 100 Myr (jMcNamara fc Nulsenl 120071 ). 

3.2. Cosmic Ray- dominated Jets 

In run JO, the jet injects into the ICM a total energy 
of Ej et — 4.21 x 10 58 erg, which is mainly kinetic. Our 
preferred jet Jl has the same amount of total jet energy, 
but it is dominated by CRs. We assume that the jet in 
run Jl has the same initial temperature Tj and velocity 
Ujet as in run JO, but has a much smaller gas density 
(rj = 10 -4 ) and consequently a very small kinetic energy. 
We choose the initial CR energy density in the jet to be 
ej cr = 2.96 x 10 -9 erg cm -3 so that £j et in run Jl is the 
same as that in run JO (see Table [T] for the jet energetics). 
To emphasize the role of CRs in the jet evolution, we 
performed another run Jl-A, which is exactly the same 
as run Jl except that the jet in this run contains no CRs 
(ej cr = 0). For reference, the central electron number 
density and temperature in the Virgo clu ster are no ~ 0.1 
cm" 3 and To ~ 1.6 keV respectively (jGhizzardi et al.l 
[2001). 

In Figures [T] and [2j we see a striking difference on the 
evolution of jets JO and Jl: Unlike the jet in run JO, 
which forms a radially-elongated cavity at large radii, 
the CR-dominated jet in run Jl forms a 'fat' low-density 
cavity near the cluster center. Due to its very low inertia 
and momentum, the light CR-dominated jet in run Jl 
decelerates much faster, and reaches a smaller distance 
than the thermal jet in run JO at a same time, as seen in 
the left panel of Figured! which also shows that the CR- 
dominated jet inflates a wider cavity which is expanded 
laterally due to the CR pressure accumulated within the 
shorter jet and cavity. As the jet turns off at tj et = 10 
Myr, the expanded jet head decelerates quickly in the 
ICM due to its low inertia, and stops moving outward 



at t ~ 20 Myr. As the jet nearly transfers all its initial 
momentum to the surrounding gas, it becomes indistin- 
guishable from the cavity, which rises slowly outward due 
to buoyancy. In contrast, the thermal jet in run JO has a 
much higher inertia and continues to move outward due 
to its initial momentum even at the end of the simulation 
(t = 100 Myr). 

At t = 70 Myr, the jet-inflated cavity in run Jl is lo- 
cated at a distance of ~ 20 kpc from the cluster center, 
having a radius of ~ 10 kpc. As the cavity rises in the 
radial direction, relatively cooler thermal gas converges 
beneath the cavity toward the jet axis, forming a ra- 
dial filament or 'thermal jet' which penetrates outward 
through the cavity. This filament is visible at t = 70 
Myr (the right-middle panel of Figure [1]) and 100 Myr 
(the middle panel of Figure [2]). This thermal jet has 
been seen in many previous simulations of X-ray cav- 
ities (e.g., IChurazov et al.l 120011 : iRevnolds etalT 120051 
Sternberg & Soker 2008b) , and i s denoted as 'cavity jet' 
byjMathews fc Brighentl (|2008at ) or 'drift' bv lPope et all 
(|201Ct ). Outward moving radial streams of cool, low- 
entropy gas are a natural post-cavity gas-dynamical de- 
velopment that probably explains cold, nearly radial fil- 
aments commonly seen in cool core clusters. 

Figure [3] shows central regions (50 x 40 kpc) of syn- 
thetic X-ray surface brightness maps for run Jl at three 
different times t = 60 (left panel), 80 (middle panel), and 
100 Myr (right panels). These X-ray maps are produced 
by integrating the radiativ e cooling rate n e niA(T, Z) 
([Sutherland fc Dopitall 19931 ) along a direction perpendic- 
ular to the jet axis (plane of the sky). Here ri[ is the total 
number density of ions and we take an average metallic- 
ity of Z = 0.4. In Figure [H X-ray cavities with radius 
~ 10 kpc are clearly seen. Interestingly, in the left and 
middle panels, the cavity is enclosed by a rim or shell 
significantly brighter than the ambient gas, a feature of- 
ten seen in real X -ray cavities (e.g.. iFabian et all 2000: 
iNulsen et al.ll2QQ2h . Dense thermal shells form naturally 
as the cavity is created, expands in the ICM and may 
become colder than the ambient gas, as suggested by ob- 
servations. 

In Figures [T] - 02 we notice that a small fraction of 
CR-filled low density gas rises up in directions perpen- 
dicular to the jet axis. As discussed in the next sub- 
section, these small cavity-like features are produced by 
backflowing cocoon gas which reaches the z = plane 
and subsequently rises in the r directions, forming small 
(toroidal) cavities perpendicular to the jet axis that have 
not been observed. In reality, the amount of backflowing 
gas reaching the z = symmetry plane may be much 
less because the jets in our simulations are initiated at 
the origin with full kpc-sized crossections, while in real- 
ity young jets are conical with much smaller crossections. 
We intend to investigate this issue in the future. 

3.3. Why do CR-dominated Jets Work? 

Why are CR-dominated jets successful in creating 
'fat' X-ray cavities near the cluster center? The CR- 
dominated jet in run Jl differs from the kinetic-energy- 
dominated jet in run JO mainly in two aspects: the dom- 
inant CR component and the lower gas density (by two 
orders of magnitude). To disentangle these two factors, 
we performed run Jl-A, where the initial jet parameters 
are the same as those in run JO, except that the thermal 
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Fig. 3. — Central regions (50 x 40 kpc) of synthetic X-ray surface brightness maps (line of sight integrated projections of the cooling rate 
n e riiA(T, Z) perpendicular to the jet axis) for run Jl at t = 60 (left panel), 80 (middle panel) and 100 Myr (right panel). The cavity is 
clearly seen as it rises in the ICM. 
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Fig. 4. — Variations of electron number density (top), the z- 
component gas velocity (middle), and pressures (bottom) along the 
r-direction (perpendicular to the jet axis) at z = 10 kpc for runs 
Jl-A (dashed; red) and Jl (other line types; black) at t = 10 Myr. 
The initial gas density and pressure profiles along the r-direction 
at t = are plotted as solid lines in the top and bottom panels, 
respectively. The dot-dashed line in the middle panel shows the 
sound speed in run Jl at at t = 10 Myr. 

gas density in the initial jet is much lower (77 = 10 -4 in- 
stead of r] = 10 -2 ), as adopted in run Jl (see Table [T] for 
the jet parameters and energetics). The bottom panels 
of Figures [T] and [2] show the jet evolution in run Jl-A, 
indicating that purely thermal jets with much lower gas 
densities decelerate much faster. In fact, thermal jets 
with lower initial gas densities have lower inertia and 
lower momentum, and are thus easier to decelerate in 
the ICM. However, unlike the CR-dominated jet in run 
Jl which creates a fat cavity resembling observed X-ray 
cavities, the under-pressured Jl-A jet creates a strongly 
elongated cavity. Furthermore, the CR-dominated cav- 
ity in run Jl is located at an even smaller cluster-centric 
distance than that in run Jl-A at a same time. Note that 
the initial jet in runs Jl and Jl-A contains the same ther- 




FlG. 5. — Variations of electron number density (top), the z- 
component gas velocity (middle), and pressures (bottom) along 
the jet axis for run Jl at t = 5 Myr (dashed) and t = 10 Myr 
(dotted). The initial gas density and pressure profiles along the 
z-axis at t = are plotted as solid lines in the top and bottom 
panels, respectively. Note that the jet is initialized in ghost zones, 
which are not plotted. 



mal component, and the difference is that only Jl con- 
tains an additional dominant CR component. Figures [T] 
and[2]show that the CR-dominated jet in run Jl creates a 
much wider cavity that experiences a larger ram pressure 
as it moves away from the cluster center. Consequently, 
the slower outward velocity of the Jl cavity quickly be- 
comes dominated by buoyancy as observations suggest, 
while the initial jet momentum may drive the continued 
outward motion of the narrow cavity in run Jl-A for a 
longer time. Thus, both the very low gas density and 
the CR pressure help the Jl jet to form a cavity near the 
cluster center. A high internal pressure from CRs in the 
Jl jet is required to create a 'fat' cavity. 

To better understand the formation of fat X-ray cav- 
ities by CR-dominated jets, we show the transverse jet 
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structure in Figure 2J which plots transverse variations 
of electron number density, the z-component gas velocity 
v z , and the gas and CR pressures along the r-direction 
(perpendicular to the jet axis) at z = 10 kpc for runs 
Jl and Jl-A at t = 10 Myr. Consider first run Jl-A 
(the dashed line in each panel). The transverse density 
profile clearly shows a bow shock at r ~ 15 kpc, which 
encloses an over-pressured cocoon of shocked cluster gas, 
a low-density cavity of shocked jet gas, and the jet mov- 
ing along the z axis. The velocity structure shown in the 
middle panel of Figure |4] indicates that the jet velocity 
drops quickly away from the jet axis, and even becomes 
negative at outer regions of the low-density cavity. Low- 
density cocoon gas with negative z-component velocities 
are backflows from the jet working surface with the am- 
bient ICM. The outer boundary of the low density cavity 
in Jl-A is located at the contact discontinuity between 
the backflowing cavity material and the shocked clus- 
ter gas, across which the gas pressure is continuous, as 
clearly seen in the bottom panel of Figure HJ By con- 
trast, for run Jl, the low-density cavity is dominated by 
CR pressure, and the total pressure (including pressures 
from both the gas and CRs) is also continuous across the 
outer boundary of the cavity. 

The low-density cavity in jet Jl is much wider than in 
Jl-A and contains thermal gas with much lower densi- 
ties, as shown in the top panel of Figure HJ Interestingly, 
the cavity gas density in run Jl (n e ~ 2 - 4 x 10 -6 cm -3 ) 
is even lower than the density of injected thermal gas at 
the jet base (n e ~ 10 -5 cm -3 for r] = 10 -4 ), which is 
a signature of the jet expansion due to the initial high- 
pressure CRs. Due to its extremely low density, the gas 
within the Jl jet and cavity will be quickly accelerated 
by any pressure gradient. In particular, the extremely 
low gas density also ensures a very high sound speed 
(eg = y (jP g + 7c P c ) 1 9) > rendering the jet velocity to be 
subsonic with respect to the jet material and the low- 
density cavity, as clearly seen in the middle panel. This 
interesting feature distinguishes very light jets from heav- 
ier jets like JO, which is supersonic with respect to both 
the jet material and its cavity. A direct consequence for 
ultra-light subsonic jets like Jl is that the total pressure 
is nearly constant throughout the cavity except near the 
origin, where the jet is quickly accelerated and expanded 
(both laterally and radially) due to the pressure gradients 
produced by the jet injection. This is shown in Figure [5j 
which plots variations of electron number density, the z- 
component gas velocity, and pressures along the jet axis 
at t = 5 Myr (dashed) and t = 10 Myr (dotted) for run 
Jl. The bottom panels of Figures |4] and [5] also show that 
the jet and cavity in run Jl are dominated by the CR 
pressure. 

The lateral jet expansion, which is essential in forming 
a fat cavity, also happens in the jet head regions near its 
working surface with the ICM. Figure [6] shows transverse 
variations of electron number density, the r-component 
gas velocity v Tl and pressures along the r-direction (per- 
pendicular to the jet axis) at z = 9.125, 8.875, and 7.875 
kpc at t = 5 Myr, when the jet at these z values is near 
the jet head, as clearly seen in the dashed lines in Figure 
[5] The middle panel of Figure [6] shows that the low- 
density jet and cavity are undergoing significant lateral 
expansion in these regions, which is driven by the small 
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Fig. 6. — Variations of electron number density (top), the r- 
component gas velocity (middle), and pressures (P g and P g + P c ) 
along the r-direction (perpendicular to the jet axis) for run Jl at 
z = 9.125 (dotted lines), 8.875 (dashed lines), and 7.875 kpc ( dot- 
dashed lines) at t = 5 Myr. The initial gas density and pressure 
profiles along the r-axis at z = 9.125 kpc are plotted as solid lines 
in the top and bottom panels, respectively. Near the jet axis, the 
gas pressure drops significantly in the low density cavity where CR 
pressure dominates. 

lateral total pressure gradient in low density gas, as seen 
in the lower panel. In fact, as the gas and CRs in the jet 
move quickly down the cluster pressure gradient, the jet 
becomes slightly over-pressured near its working surface, 
producing the lateral pressure gradient. The resulting 
cavity is fatter when the injected jet pressure is larger at 
the jet base. 

For a direct comparison with the kinetic-energy domi- 
nated JO jet, our initial parameters for the CR-dominated 
Jl jet are chosen so that it has the same energy flux 
as JO, but mainly in the form of CRs. To this end, 
we chose the initial CR energy density in the jet to be 
ej cr = 2.96 x 10 -9 erg cm -3 , which makes the jet over- 
pressured at the jet base (the ratio of the jet-to-ICM 
pressure is about 2). In principle, the initial jet pres- 
sure is determined by physical processes directly related 
to the central super massive black hole, and is unlikely 
to equal to the pressure of central hot ICM. The jet may 
become over-pressured as it quickly moves down the clus- 
ter pressure gradient. Moreover, CRs may be produced 
within the jet, rendering the jet suddenly over-pressured 
(e.g., Laing et al. 2006). Run Jl thus represents these 
cases. 

We also consider another run J2 in which the jet is ini- 
tially in pressure equilibrium with the ambient ICM at 
its base. The J2 jet differs from Jl only in having a lower 
CR energy density at the jet base, ej cr = 1.52 x 10 -9 erg 
cm -3 . The total jet energy flux in run J2 is thus less 
than that in run Jl (see Table 1). Figure [71 shows a cen- 
tral 2D slice (50 x 20 kpc) of electron number density in 
logarithmic scale in run J2 at time t = 60 Myr. As in run 
Jl, the CR-dominated jet in run J2 produces a fat X-ray 
cavity near the cluster center. The jet and cavity struc- 
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ture is very similar to that in run Jl. In contrast, the 
jet in run Jl-A without CRs produces a strongly radially 
elongated cavity, suggesting that a minimum amount of 
CR component is required to produce a fat X-ray cavity 
near the cluster center. 

Successful CR-dominated jets discussed above have ex- 
tremely low initial densities (77 = 10 -4 ) and the jet en- 
ergy flux and pressure are mainly carried by CRs. We 
have also performed a few simulations of CR-dominated 
jets with larger gas densities 77 = 10 -2 (i.e., with jet elec- 
tron number densities rjno ~ 10 -3 cm -3 ). By increasing 
the CR pressure within the jet, we find that as expected, 
the cavity produced rises in the ICM much slower than 
in run JO, where no CRs are introduced. However, due to 
the higher jet inertia and momentum, the cavity formed 
is usually significantly elongated in the jet direction, sim- 
ilar to run JO. This can be seen in Figure [8j which illus- 
trates a central 2D slice (50 x 20 kpc) of log (n e /cm -3 ) 
in a typical run J3 at time t = 70 Myr (see Table 1 for 
simulation parameters in this run). Although the back- 
flowed gas forms a fat cavity that rises slowly in the ICM 
(mainly due to buoyancy), a large fraction of the jet ma- 
terial, in particular the jet head, still moves quite fast, 
driven by the large initial jet momentum. The cavity 
formed is significantly radially elongated, strikingly dif- 
ferent from the cavity produced in run Jl (Figures Q] and 
[2]). Thus our preliminary calculations suggest that ob- 
served X-ray cavities are likely to have been produced by 
CR-dominated light jets with 77 <C 0.01. Jet simulations 
involve a large parameter space, and we leave a more 
thorough parameter study to future work, which we an- 
ticipate will also be possible to explain observed X-ray 
cavities and radio bubbles spannin g; a large range of size, 
clust er-centric distance, and age (McNa mara fc Nulsenl 
l2007h . 

The main dynamical effect of CRs is their contribution 
to the ICM pressure while introducing negligible inertia. 
The same effect may also be achieved with a purely ther- 
mal jet by proportionally increasing the gas temperature 
while extremely reducing the gas density. However, to 
reach a jet-to-ICM density ratio of 77 = 10 -4 , the jet tem- 
perature should be at least - 10 4 keV (i.e. T > 10 11 K), 
and gas with such high temperatures is already mildly 
relativistic [electron gamma- factor > 30(/cbT/10 4 keV)]. 
Radio synchrotron emission from FR I jets and some X- 
ray cavities indicates a significant population of CR elec- 
trons with higher energies (gamma-factor of ~ 10 3 — 10 5 ). 
In our model, the dynamics of CRs is described by e c , 
the integrated energy density over the CR energy distri- 
bution. Our model is applicable to relativistic CRs with 
any energy spectrum or any content (i.e., any combina- 
tion of CR electrons and protons). A proper consider- 
ation of these details is necessary to calculate the radio 
emission from CRs and this is an important problem for 
the future. 

4. CONCLUSION AND IMPLICATION 

AGN feedback has been frequently invoked to explain 
the cut-off of the galaxy luminosity function at the bright 
end and the suppression of cooling flows in giant ellip- 
tical galaxies, galaxy groups, and clusters. Many X-ray 
deficient cavities detected by high-resolution X-ray ob- 
servations provide strong evidence for AGN feedback in 
galaxy clusters. Most cavities are roughly spherical or 
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Fig. 7. — Slice of log (n e /cm~ 3 ) in run J2 at time t = 60 Myr. 
Contours show CR energy density e cr in units of 10 -10 erg cm -3 
at four levels 0.01, 0.1, 0.5, 1.0. 
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Fig. 8. — Slice of log (n e /cm -3 ) in run J3 at time t = 70 Myr. 
Contours show CR energy density e cr in units of 10 -10 erg cm -3 
at four levels 0.01, 0.1, 0.5, 1.0. 

slightly elongated in the tangential direction, and those 
detected are commonly located near cluster centers (off- 
set by a few to tens kpc), having spatial sizes of a few 
to tens kpc. However, in numerical simulations it has 
been difficult to reproduce such fat cavities. Hot purely 
thermal jets frequently adopted in previous simulations 
usually penetrate deeply and rapidly through the ICM, 
forming radially-elongated cavities at much larger radii. 

In this paper, for the first time, we present numerical 
simulations of CR-dominated AGN jets and follow their 
dynamical evolution as they form cavities in the ICM. Us- 
ing two-dimensional axisymmetric calculations, we show 
that CR-dominated jets naturally create fat cavities near 
cluster centers, resembling observed cavities very well. 
The key feature of successful CR-dominated jets in our 
model is that they are also extremely light - the jet- 
to-ICM density ratio is 10 -4 in our main run Jl. The 
low inertia and momentum in these jets ensure that they 
are quickly decelerated in the ICM. The energy fluxes 
of these jets are dominated by CR energy, and the CR 
pressure drives the lateral jet expansion, which assists 
the jet deceleration and is essential in forming fat lobes. 
The resulting lobes are mainly formed by jet backflows, 
and supported by the CR pressure, which also displaces 
ambient thermal gas, naturally creating X-ray deficient 
cavities filled with synchrotron emitting CR electrons. 
The initial low jet momentum is quickly transferred to 
the ambient post-shock cluster gas, and the X-ray cavi- 
ties then rise buoyantly in the ICM. 

Our scenario in which X-ray cavities evolve from CR- 
dominated AGN jets is supported by observations of CR- 
filled X-ray cavities in elliptical galaxies, galaxy groups 
and clusters which are often associated with synchrotron- 
emitting radio jets and spatially coi ncident with radio 
bubbles (|McNamara fe Nulsenl I2QQ7D . High-resolution 
multi- frequency radio observations of FR I jets, which are 
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probably responsible for forming X-ray cavities, indicate 
that these jets decelerate rapidly and produce strong syn- 
chrotron emission (flaring) at typical distances of a few 
kpc from central nuclei (e.g., Laing et "aTll2QQ6f ). These 
observations thus also point toward a significant CR elec- 
tron population in FR I jets on kpc scales. 

Our calculations of cavity formation imply that, to cre- 
ate observed fat X-ray cavities near cluster centers, AGN 
(FR I) jets must be (1) very light, containing a very small 
thermal component, and (2) dominated by CR pressure. 
We thus derive new constraints on the composition of FR 
I jets solely from their ability to create fat X-ray cavities 
(radio bubbles) as observed in elliptical galaxies, groups 



and clusters of galaxies. The thermodynamical impact 
of CR-dominated jets on the ICM, e.g., their efficiency 
in heating cooling flows and spreading metals, may be 
different from that of purely thermal jets, and this will 
be investigated in future studies. 



We thank the anonymous referee for a positive and 
detailed report, which was very informative and helpful. 
Studies of CR-dominated AGN jets and the evolution of 
hot cluster gas at UC Santa Cruz are supported by NSF 
and NASA grants for which we are very grateful. 
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